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Recently, the effects of spin-orbit coupling (SOC) in correlated materials have become 
one of the most actively studied subjects in condensed matter physics, as correlations and 
SOC together can lead to the discovery of new phases. Among candidate materials, iridium 
oxides (iridates) have been an excellent playground to uncover such novel phenomena. In 
this review, we discuss recent progress in iridates and related materials, focusing on the 
basic concepts, relevant microscopic Hamiltonians, and unusual properties of iridates in 
perovskite- and honeycomb-based structures. Perspectives on SOC and correlation physics 
beyond iridates are also discussed. 
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I. INTRODUCTION 

Spin-orbit coupling (SOC) is a relativistic effect that links the orbital and spin angular momenta 
of an electron. Though suppressed by the fine structure constant, the electric fields near the nuclei 
of atoms with a large number of protons can render this interaction significant. The natural place 
to find significant SOC is thus in atoms with high atomic numbers, moving down the rows of the 
periodic table into the heavier elements. 

The effects of SOC in materials with such heavy atoms have been studied intensively in the 
context of semi-conductors. In these weakly correlated materials, SOC entangles the crystal mo¬ 
mentum and spin of the electron, locking the kinetic and internal degrees of freedom together. 
This leads to a number of intriguing phenomena, particularly in transport; examples include the 
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anomalous Hall effect (AHE) and the control of spin currents being applied in the field of spin- 
tronics. More recently it has been found that SOC plays an essential role in the fast growing field 
of topological insulators (TI) and metals. Many, if not all, of the experimental examples of such 
topological phases have been found in these types of heavy semi-conductors. 

While the effects of SOC in weakly correlated materials described above has been thoroughly 
considered, its importance in more strongly correlated transition metal materials remains less de¬ 
veloped. There are several energy scales to consider in such materials: the atomic interactions 
(schematically) on-site Hubbard interaction U, Hund’s coupling J H , the SOC A, the crystal field 
A and the electron kinetic energy described by hopping integral t. Much theoretical and experi¬ 
mental effort has been brought to bear on 3d transition metals such as high temperature cuprates, 
manganites, and vanadium oxides. In these compounds the atomic interactions, crystal field and 
kinetic terms can generally compete, though the SOC remains small. In heavy transition metals 
such as those with 5 d and even Ad electrons SOC is significant as well, and so all of these energy 
scales can be comparable. As one moves through the different heavy transition metal materials, 
small changes in these details can tip the scales, revealing a surprisingly rich family of behaviours. 
This can be contrasted with 4/ or 5/ electrons in the lanthanides or actinides, where the electron 
interactions are dominant, followed by the SOC and then the crystal fields. 

While intensive studies have recently been undertaken, we are still far from complete under¬ 
standing on the combined effects of SOC and electronic correlation. Material examples of Ad and 
5 d systems can be found for a variety of two- and three-dimensional lattices with varying degrees 
of frustration, covering the range from weakly to strongly correlated. An excellent overview of this 
physics is provided in a recent review 1 , with a focus on examples from the pyrochlore iridates and 
the 5 d double-perovskite magnets. In this review we will focus on magnetism, topological phases 
and superconductivity induced by the competition between strong SOC and electronic correlation. 
We illustrate these concepts with examples drawn from the perovskite and honeycomb iridium 
oxides. After a basic introduction to the relevant models in Sec. II, we begin with perovskite 
iridates in Sec. III. Drawing analogies to the related cuprates and ruthenates, we discuss proposals 
for realizing topological and superconducting phases in the materials. In Sec. IV we tackle the 
burgeoning field of Kitaev magnetism. In particular, we offer a summary of the current status of 
two- and three-dimensional honeycomb iridates. Concluding remarks and outlook are provided in 
Sec. V. 


II. BACKGROUND 

The building blocks of our discussion are the atomic states of the partially filled Ad- or 5 d ions. 
In solids of interest, these states are split by a predominantly octahedral crystal field potential into 
a t 2g triplet and e g doublet, as shown in Fig. 1. The energy gap to the e g doublet is large, so 
these states can be safely ignored when we consider electron filling less than six. We are primarily 
interested in the d 5 configuration, which can be regarded as single hole in one of the t 2g states. 
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When projected into this manifold, the angular momentum of the the d electrons is mapped to a 
set of effective / = 1 angular momentum operators, -L. The large SOC then acts within the t 2g 
manifold as -AL ■ S where L is an effective / = 1 angular momentum and S is the spin. Using 
the rules of addition of angular momenta we see that SOC splits the t 2g multiplet into an effective 
j - 1/2 doublet and effective j = 3/2 quartet as show in Fig. 1. The j = 3/2 states are lower in 
energy and separated from the j = 1/2 states by a gap of 3A/2. Written in terms of the t 2g states 
one has 
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where 1 yz, ±), |jcz, ±), | xy, ±) are the t 2g states and spin up and down correspond to ±. The im¬ 
portant role played by SOC can be seen in the entanglement of the spin and orbital degrees of 
freedom in these wave-functions, as illustrated in Fig. 1. The states of this j = 1/2 doublet and 
its pseudo-spin operators J will be a common thread in our discussion of the family of iridium 
oxides and ruthenium materials. While strongly spin-orbit entangled, it remains isotropic, with a 
^-factor of -2. Explicitly, the magnetic moment operator p = p B (L + 2S) becomes -2p B J when 
projected into the j = 1/2 states. The anisotropy of these systems thus manifests predominantly 
in the coupling of j = 1/2 moments , not in the single-ion properties. 

The interactions of j = 1 /2 electrons are determined by the atomic interactions of the free ion 
projected into the t 2g manifold and the kinetic terms, or hoppings of the t 2g electrons. Schematically 
we can write the general multi-orbital model 


ZZZ^(4^+ h - c )+E 


i j ap cr 


U-3J , 


H 


{Ni - 5 ) 2 - 2J h Si 




H l2- AL, • S; 


( 2 ) 


where Ni is the total number operator, S, is the total spin operator and L, is the total pseudo- 
angular momentum operator at site r ( . The d] a(T operator creates a t 2g electron in orbital a = yz, xz 
or xy with spin cr j, at site r, and //f are the tight-binding hopping parameters. The hopping 
amplitudes f/:‘ will have contributions coming from direct d-d overlap, as well as from processes 
through the intermediate oxygen atoms. As these kinetic terms are strongly material dependent 
in both structure and scale, we save discussion of their details for the sections devoted to each 
specific class of material. 

The local atomic interactions are more generic and can be expressed in terms of two parameters ' : 
the Coulomb repulsion U and the Hund’s coupling J H . We have added a chemical potential to 
favour the N = 5 state relevant for Ir 4+ and Ru 3+ and have fixed the inter-orbital repulsion V to the 
free ion value U-2J H . Typically one expects U ~ 2 eV , J H ~ 0.2 eV and 3A/2 ~ 0.4-0.5 eV for 
an Ir 4+ ion. On the other hand for a Ru 3+ ion, U and J H are stronger while atomic SOC is weaker 
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FIG. 1: (a) Octahedral geometry of transition metal site illustrated for an iridium oxide, (b) 
Splitting of the 4 d or 5 d levels by octahedral crystal fields A 0h and by SOC A into j = 1/2 and 
j = 3/2 levels, (c) Illustration of the atomic j = 1/2 wave-functions. The composition of the 
j = 1/2 states with the spin-f charge density is shown in red and the spin-| charge density in blue. 

than the Ir 4+ case. Due to screening effects in a solid, these parameters will be renormalized, and 
the free ion relation U' - U - 2 J H can be violated. Typically, one expects the Coulomb integral U 
to be more strongly screened than the Hund’s coupling J H ' . 

There are several limiting regimes where the physics of this model is particularly clear. In Sec. 
Ill we will consider the limit where the interactions between the j - 1/2 and j - 3/2 electrons can 
be neglected. This yields an effective single-band Hubbard model for the j = 1/2 electrons. With 
sufficiently weak correlations, this approximation can be useful as a starting point for itinerant 
single band j = 1/2 systems. In the strongly correlated regime, the j = 1/2 electrons localize 
and one is left with a pseudo-spin model. By projecting into only the j = 1/2 states one loses 
the j = 3/2 excited states and the effects of Hund’s coupling. As we will see in Sec. IV, for the 
honeycomb iridates, this leading term may cancel , so the contribution of the virtual processes that 
go through j = 3/2 states become the largest interactions. 


III. MAGNETISM, TOPOLOGICAL PHASES AND SUPERCONDUCTIVITY IN 

PEROVSKITE IRIDATES 

Historically, the field began with the study of Sr 2 Ir0 4 and Sr 2 Ru0 4 as isostructural analogues 
of La 2 Cu0 4 , the parent compound of the cuprate superconductors. While Sr 2 Ir0 4 was first syn¬ 
thesized in 1957 1 , early work focused mainly on chemical aspects, with variety of ternary and 
quaternary oxides of iridium and ruthenium reported in the following decades. 1 More serious 
interest in Sr 2 Ir0 4 was piqued when superconductivity was discovered in Sr 2 Ru0 4 \ but an insu- 
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FIG. 2: Crystal structures of the perovskite iridates from the Ruddlesden-Popper series 
Sr n+ ilr n 0 3n+ i. These include the quasi-two dimensional (a) single-layer Sr 2 Ir0 4 . and (b) 
bilayer Sr 3 Ir 2 0 7 and the three dimensional limit (c) orthorhombic Srlr0 3 • . 

lating state with small ferromagnetic moment was found instead - . Interest was renewed by the 
discovery that the interplay between SOC and electronic correlations were driving the physics in 
this material - . 

In this section, our primary focus will be on magnetism in Sr 2 Ir0 4 and Sr 3 Ir 2 0 7 and the pos¬ 
sible realization of novel topological phases and superconductivity. We start with a review of the 
experimental and theoretical progress in understanding perovskite iridates such as the Ruddlesden- 
Popper series of Sr n+ ilr n 0 3n+ i through a discussion of the common building block of corner-shared 
octahedra in both the itinerant and localized limits. Next, we discuss some candidate topological 
phases built from iridium oxide heterostructures. Finally, we discuss the effects of charge doping, 
with an eye toward the possibility of superconductivity in Sr 2 Ir0 4 . 


A. Magnetism 

In the single and bilayer perovskites Sr 2 Ir0 4 and Sr 3 Ir 2 0 7 6 evidence from the optical gap 
to the magnetic ordering transition - ' have indicated that both are Mott insulators, though 
Sr 3 Ir 2 0 7 being substantially weaker than Sr 2 Ir0 4 . The modest on-site Hubbard interaction U 
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proves sufficient to localize the electrons due to the narrow, spin-orbit coupled j = 1/2 bands. In 
Sr 2 lr 04 the magnetic ordering is an in-plane, canted antiferromagnet • following the staggered 
rotation of the oxygen octahedra, as illustrated in Fig. 3a. Though the net ferromagnetic moment 
of this state was observed in early studies* ’ ,14 , the crucial role played by strong SOC was only 
understood fairly recently . The moment directions in Sr 3 Ir 2 07 are very different, forming a 
collinear arrangement perpendicular to the iridium planes^ as shown in Fig. 3c. 

The bulk trilayer and higher compounds {n > 3) are unstable under ambient pressure , includ¬ 
ing the three-dimensional limit Srlr0 3 which takes a post-perovskite structure . An orthorhombic 
perovskite Srlr0 3 can be stabilized via the application of pressure; this is a semi-metal and shows 
no evidence of magnetic ordering ■ . One can then posit a metal-insulator transition (MIT) 
exists as a function of layer count near n c ~ 3, placing Sr 2 Ir0 4 (n = 1) and Sr 3 Ir 2 0 7 (n = 2) on the 
insulating side and orthorhombic Srlr0 3 (n = oo) on the metallic side. Due to the proximity to this 
putative MIT there has been some debate on the applicability of the strongly coupled Mott picture 
as compared to a more weakly coupled Slater picture ~ . Since the weak and strong-coupling 
regimes give the same magnetic orderings , 2 this distinction may be somewhat academic at low 
temperature. Near the transition the Mott and Slater regimes could be distinguishable through the 
behaviour of the charge gap. 

Key questions to address in these compounds concern both the magnetic ordering and the MIT. 
In particular, we wish to understand what sets the direction of the ordered moments in Sr 2 Ir0 4 , and 
Sr 3 Ir 2 0 7 . Beyond this we would like to establish a framework understanding both the itinerant and 
localized limits of this family of compounds. This will provide useful guidance for our discussion 
of engineered phases constructed from the same basic building block in Sec. HID. 


B. Corner shared octahedra with strong spin-orbit coupling 

The perovskite iridates are built from corner-shared octahedra as shown for Sr 2 Ir0 4 , Sr 3 Ir 2 0 7 
and orthorhombic Srlr0 3 in Fig. 2. The observed MIT as a function of layer count makes 
it necessary to discuss both the weak and strong correlation limits. Here, we present a brief 
account of the j - 1/2 spin models at strong coupling, as well as the tight-binding models at weak 
coupling that is applicable to all networks of corner-shared octahedra. While the real materials 
have structural distortions, we begin with the idealized, undistorted limit. 

One simple approach to attacking the model in Eq. (2) is to assume the j = 1/2 and j = 3/2 
bands to be well separated, neglecting the inter-band interactions. This projects the multi-orbital 
problem of Eq. (2) into a single-band model of j = 1/2 electrons with an effective hopping and 
Hubbard interaction 

2 2 t{ j +c 5« c i «) +Udr 2 ni+n ‘-’ (3) 

ij a=± i 

where c\ + creates a j z = +1/2 electron at site r, and t/ eff ~ ( U+2U')/3. Inversion and time-reversal 
symmetry forbid pseudo-spin dependent hoppings. This type of itinerant model is a useful starting 
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point to describe the more three-dimensional compounds, such as those based on Srlr 03 . Due to 
the increased dimensionality, the bandwidth of the j = 1/2 states is larger than the quasi-two- 
dimensional compounds and becomes comparable to the effective Coulomb C/ eff . 

The magnetic ordering in these models is most easily approached from the strong Mott limit 
where U e s » t. Assuming the nearest-neighbour hoppings are dominant, this yields an antiferro¬ 
magnetic Heisenberg model with an effective SU(2) pseudo-spin rotation symmetry 

— Yj -J 


(4) 


where J, is the j = 1/2 pseudo-spin at site r, and t is the nearest-neighbour j = 1/2 hopping. 
Going back to Eq. (2) and including the j = 3/2 excited states or the Hund’s coupling spoils 
this accidental symmetry and generates all symmetry allowed terms. For an ideal corner-shared, 
180° bond these include a four-fold rotation about the bond direction, and several reflections. The 
symmetry allowed terms include a Heisenberg exchange J and a compass-like, pseudo-dipolar 
coupling K 

J] [7J, ■ Jj + K(rij ■ J;)(r, 7 • J ; ), ] (5) 

Hi) 

where r, ; = r y --r,- is the bond direction. In concrete models- one finds K/J oc J H /U, so we expect 
the Heisenberg term to dominant, with only a small pseudo-dipolar coupling K J. 

The structure of this simple model provides a partial answer to the experimental questions 
introduced in Sec. Ill A. Namely, the approximate SU(2) symmetry renders the final pinning 
of the magnetic moments sensitive to SU(2) breaking perturbations, such as the pseudo-dipolar 
coupling K. Due to this sensitivity, to understand the fate of these moments we must understand 
the complete set of structural distortions present in the perovskite iridates. 


C. Distortions and octahedral rotations 

Two relevant types of distortions in the perovskite iridates are tetragonal distortion and octahe¬ 
dral rotations. The former arises in the layered compounds such as Sr 2 Ir 04 and Sr 3 Ir 2 07 , while 
octahedral rotations are present in the entire Ruddlesden-Popper series. Due to the layered struc¬ 
ture, the octahedra in Sr 2 Ir0 4 and Sr 3 Ir 2 0 7 are elongated along the inter-layer direction, [001] l_: ' . 
Such local distortions can be encapsulated in a residual crystal field potential that split the j - 3/2 
levels into two Kramers doublets which remain well separated from the higher-lying doublet for 
small distortion . This highest lying doublet plays the role the j = 1/2 does in the ideal case and 
defines our pseudo-spin. In their simplest form, these distortions manifest as a compression or 
elongation along some direction n of the octahedron, modeled as 

Vfi = A(n • L) 2 , (6) 

where L is the orbital angular momentum projected into the t 2g levels and A > 0 corresponds to 
compression and A < 0 to elongation. For tetragonal distortion n is along one of the cubic axes, 
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FIG. 3: Illustration of octahedral rotations and magnetic canted antiferromagnetic order in 
Sr 2 Ir0 4 • . (a) The octahedral rotations and moment directions in a single Ir0 2 plane of Sr 2 Ir0 4 . 
(b-c) Illustration of the stacking of planes and the ordering in (b) Sr 2 Ir0 4 ’ and (c) Sr 3 Ir 2 0 7 . 


x, y or z. Including the SOC, the atomic energy levels can be straightforwardly found, with the 
highest lying doublet being a mixture of the j = 1/2 and j = 3/2 states as 


|+) = cos 8 



+ sin# 



(7) 


where the mixing angle is tan(2#) = 4 V2A/(2A + 9A). Including these distortions removes the 
effective four-fold rotation symmetry, allowing an additional Ising-like anisotropy ~ r zz J:Jj on 
each bond in the plane perpendicular to the distortion . 

In addition, octahedral rotations are present and break some of the rotational and translational 
symmetry of the ideal lattice, enlarging the unit cell. In Sr 2 Ir0 4 and Sr 3 Ir 2 0 7 the octahedra are 
rotated about the [001] axis in a staggered fashion 1-8 , characterized by a deviation of the bond 
angle of ~ 11° = do from the ideal 180°, as illustrated in Fig. 3a. We leave the more involved 
octahedral tilting of orthorhombic Srlr0 3 to the literature ■ . Such rotations lower the bond sym¬ 
metry, breaking the inversion symmetry about the bond center and allow pseudo-spin dependent 
hoppings and Dzyaloshinskii-Moriya (DM) exchange ■ . Schematically, a tight-binding model 
for the j = 1/2 electrons in a single Ir0 2 plane that includes such effects is 


Z 

Hi) 


Zfe“ +hx -) + ^Z (_1) '( 


cr,c 


z Cj a ' 


+ 


h.c.) 


( 8 ) 


where (-1)' = ±1 is a staggered sign, as illustrated in Fig. 3a and one expects t ~ cos 2 0 o and 
t z ~ sin 20 o - The pseudo-spin dependent term t z can be removed by rotating the pseudo-spin at 
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each site by an angle (-1 ) ! 0 about the z axis, where tan 20 = t z /t 2H . In these rotated axes the 
effective spin model is once again purely Heisenberg 2 , with J ~ 4(t 2 + t 2 )/U c a, without the DM 
or symmetric anisotropic exchanges that were expected in terms of the original pseudo-spins. 

This is an accident of the strict j = 1/2 limit ; the symmetric anisotropies, such as T z , intro¬ 
duced by tetragonal distortion, the pseudo-dipolar K or the inter-layer anisotropies in the bilayer 
case will not be removed. These break the pseudo-SU(2) symmetry and thus set a preferred di¬ 
rection or set of directions in pseudo-spin space. An appropriate choice of these terms can 
pin the moment in a direction that is staggered following the octahedral rotations as is found in 
Sr 2 Ir0 4 . The case of Sr 3 Ir 2 07 is less clear, though the moment orientation being very different is 
not surprising in light of the above discussion. However, it is likely this compound lies far from 
pseudo-SU(2) invariant limit given the very large gap seen in the spin-wave spectrum . While 
consensus has not yet been reached, several proposals have put forth mechanisms to explain this 
collinear order by including large anisotropic couplings ■ . 


D. Topological phases 

While in the above sections we started from the Mott insulating limit, in moving to orthorhom¬ 
bic Srlr0 3 we will start from the itinerant limit. As bulk Srlr0 3 is a semi-metal with large SOC, 
variants of this compound may be a good place to look for topological insulators or metals. A 
band insulator is deemed a topological insulator (TI), if it cannot be smoothly deformed into a 
decoupled, atomic insulator in such a way that preserves time-reversal symmetry 1,32 . In two or 
three dimensions this classification into trivial and non-trivial is complete without any further 
subdivisions ’ . TIs have a number of very interesting properties such as gapless edge states and 
unusual magneto-electric response • . As discussed in the introduction, most examples of TIs are 
found in weakly correlated systems. Some proposals have been put forward to look for TIs in tran¬ 
sition metal oxides (TMOs): these promise to not only be more robust, but provide a playground 
to explore the effects of interactions on such topological phases 1,3 '. 

Among the iridates, it was suggested that Na 2 Ir0 3 could realize a TI with the pseudo-spin 
j - 1/2 band acquiring a non-trivial Z 2 topological index' . While there were some early attempts 
at explaining the magnetism of Na 2 Ir0 3 within a band-type picture , experimentally it appears 
that Na 2 Ir0 3 is a Mott insulator. In three dimensions, the pyrochlore iridates were also put forth 
as candidates to realize a strong TI . Another promising approach is to design the materials of 
interest by adopting techniques developed in studies of oxide heterostructures 36 . These methods 
allow for a wide variety of lattice geometries with a greater control over structural distortions and 
impurity content than is available in bulk samples. Heterostructures built from 4 d or 5 d TMOs, 
with their large intrinsic SOC, are thus likely to offer new and interesting phases for study. 

A simple example of such physics can be illustrated via bilayers of Srlr0 3 grown along the 
[111] direction ■ . In these bilayers the Ir atoms form a buckled honeycomb lattice of corner- 
shared octahedra, as shown in Fig 4a. As in the proposals to realize a TI in Na 2 Ir0 3 , we have an 
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FIG. 4: (a) An ideal bilayer of Srlr 03 showing the honeycomb lattice formed by the Ir ions. First 
nearest neighbours are shown by solid lines while second nearest-neighbour bonds are dashed. 
The arrow indicates the hopping direction, while the colour indicates whether d= x, y or z. (b) 
An ideal superlattice [(SrIr0 3 )3,SrTi0 3 ] along the [001] direction. 

itinerant j = 1/2 model of the form 

- ? Z Z ( c L c * - hx -)+ f Z Z K K - 6 ■ *") +h - c l +e Z ( - |) '^- (9) 

<«/> « «' 7 » a P ia 

where cj a is a creation operator of j = 1/2 electron with the pseudo-spin a = j z = +1/2 at 
site r ( . The sums over (..) and ((..)) denote the nearest neighbour (NN) and next nearest neighbour 
(NNN) bonds, respectively. On the NNN bonds the absence of inversion symmetry allows pseudo¬ 
spin dependent terms to appear parametrized by d,- ; = x, y or z as illustrated in Fig. 4a. Due to 
the buckling of the honeycomb lattice, there is an atomic potential e staggered between the two 
sublattices. With t' = e = 0, this system is a semi-metal with Dirac cones, as in graphene, while 
finite f or e gaps out these cones giving a band insulator. When the atomic potential is sufficiently 
small, this model supports a TI with a non-trivial Z 2 index and the associated gapless edge modes. 

The simple example presented above only begins to scratch the surface of what is possible in 
these types of engineered systems. One particular class of system that has been systematically 
studied is a series of artificial superlattices [(Srlr0 3 ) m , SrTi0 3 ] where m is integer, grown along 
the [001]-axis atop a substrate SrTi0 3 . These superlattices consist of atomically thin slices of 
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Srlr 03 separated by layers of insulating SrTi 03 . While bulk Srlr 03 forms a post-perovskite struc¬ 
ture under ambient conditions 9 , when prepared in thin films or superlattices 23 , the SrTi 03 layers 
stabilize the orthorhombic perovskite structure. By tuning the number of Ir0 2 layers, an insulator 
was achieved for m = 1 single and m = 2 bilayer superlattices, and trilayer m = 3 illustrated in 
Fig. 4b sits at the verge of MIT driven by the layer number m . While these m = 1 and in = 2 
superlattices are topologically trivial j = 1/2 insulators, a theoretical study on single and bilayer 
Srlr0 3 interleaved with a band insulator with orthorhombic perovskite structure (such as GcSc0 3 , 
CaTi0 3 , SrZr0 3 or SrHf0 3 ) suggests a rich phase diagram including topological magnetic insula¬ 
tors, topological crystalline insulators and topological valley insulators 40 . 

The stability of some topologically protected states hinges the presence of symmetries. A richer 
topological classification has been recently uncovered by considering the phases protected by the 
spatial symmetries of the crystal. These topological crystalline insulators share many of the same 
features as Z 2 TIs, such as protected gapless surface states 41 . In contrast, in the prototypical Weyl 
semi-metal, the gapless surface states remain protected even in the absence of any symmetry. 
Given the finer classification that exists for insulators, one can ask if there are topological met¬ 
als that are protected by spatial symmetries. A large class of such topological metals have since 
been classified 3 ; these go beyond the Weyl semi-metal and have surface states protected by ei¬ 
ther global symmetries or crystal lattice symmetries, or combination of the two. Depending on 
the symmetry properties and the dimension of the Fermi surface, these surface states can form 
Dirac cones, flat bands or Fermi arc states. . Among the iridates, the three dimensional per¬ 
ovskite iridates Alr0 3 with A an alkaline earth metal have been recently proposed as an example 
of topological metals protected by spatial symmetries. These topological crystalline metals are 
analogous to the topological crystalline insulators, with crystal symmetry responsible for the pro¬ 
tected surface states. While the bulk states show a nodal of ring of gapless excitations, protected 
by time-reversal and spatial symmetries, the associated two-dimensional surface states are flat in 
one direction, while linearly dispersing in the other. If one adds external symmetry breaking terms 
such as those that break time-reversal, mirror, or glide symmetry, this topological metal acts as a 
seed for a rich family of topological phases, such as a Weyl semi-metal, strong or weak topological 
insulator, or topological magnetic insulator 44 . 


E. Superconductivity 

We now return to the original motivation for the study of the iridates and Sr 2 Ir 04 in particular: 
possible links to the physics of high-temperature superconductivity and the realization of exotic 
superconducting pairing symmetry in doped spin-orbit Mott insulators. In multi-orbital systems 
it is a challenge to understand the role played by the orbital degrees of freedom in the microscopic 
mechanism of superconductivity. The superconductor Sr 2 Ru 04 provides a relevant example, with 
all three t 2g orbitals important in constructing Fermi surface. Even such a 4 d transition metal 
SCs with intermediate SOC, multi-orbital interactions such as Hunds coupling can be important. 
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FIG. 5: (a) An illustration of the energy levels for the iridium and oxygen atoms. A hole on the 
oxygen is more costly than one in the j = 3/2 states, (b) Schematic phase diagram for Sr 2 Ir0 4 as 
a function of filling n and Hund’s coupling J H , based on Ref. [47]. A variety of phases appear 
such as an antiferromagnetic insulator (AF-I) and metal (AF-M), a paramagnetic metal (PM-M), 
a ferromagnetic metal (FM-M) as well as possibly two distinct superconducting (SC) phases. 


Moving to 5 d transition metals, the large SOC can no longer be ignored and the system may 
display a complex combination of spin-singlet and spin-triplet SC order parameters. As the SOC 
and other electronic interactions such as Hund’s coupling become comparable, the determination 
of the ground states in such multi-orbital systems is highly non-trivial. In particular, we will 
discuss how the physics of doping in Sr 2 Ir 04 is both similar to, but different from the case of 
cuprate superconductors. We then discuss some recent theoretical work that explore possible 
pairing symmetries and mechanisms for superconductors arising from doped Mott insulators with 
strong SOC. From there, we review the on-going experimental search for superconductivity in 
doped iridates, focusing in particular on progress in doping Sr 2 Ir0 4 . 

The physics of doping electrons into these j = 1/2 Mott insulators is familiar. Given the filled 
2 p orbitals of the neighbouring oxygens and the filled j = 3/2 states, the least costly place to put 
an extra electron is in the j = 1/2 states themselves. The penalty is the Coulomb energy U ; much 
smaller than filling any of the higher-lying electronic states. We caution that this is purely at the 
atomic level and electronic structure effects can alter this identification. For example for Sr 2 Ir0 4 
it has been argued that the electron doping is similar to the hole-doped cuprates due to an opposite 
sign in the next-nearest neighbour hopping integral . However when iridates are hole-doped, the 
multi-orbital nature can become active and it may not be this simple. The key distinction is in 
the energy gap to the nearest 5 d states, the j = 3/2 levels. In the cuprates, tetragonal distortion 
separates the ground state doublet from the nearest 3d levels by a gap of ~ 1 - 2 eV . This is larger 
or comparable to the cost of putting the hole onto the neighbouring oxygens once hybridization 
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effects are taken into account. In the case of the iridium oxides the j = 3/2 states are only ~ 0.5 eV 
away while the oxygens are ~ 1 - 2 eV. It is then favourable to put the hole not on the oxygens, 
but on Ir 4+ itself in the j = 3/2 states. Hund’s coupling between the j = 1/2 and j = 3/2 states 
can further stabili z e these holes. This is schematically summarized in Fig. 5a. A single-orbital 
j - 1/2 picture, enabled by the formation of a Zhang-Rice singlet, may no longer be tenable and 
one must move beyond a t-J type model and include aspects of the full multi-orbital problem. This 
can be seen clearly in how these same j = 3/2 excited states appear in the localized spin model for 
the half-filled d 5 case, where they generate bond-dependent exchanges as discussed in Sec. TUB 
and will be essential in Sec. IV. 

The classification of pairing terms is also modified when SOC is strong; as spin is no longer a 
conserved quantity spin-singlet or spin-triplet lose their meaning. However a version of this can 
be recovered if the strong SOC leads to a single FS with j = 1/2 character and the multi-orbital 
problem can be reduced to an effective single band model. If the pseudo-spin dependent terms can 
be made small, one can define effective pseudo-spin-singlet and pseudo-spin-triplet SC pairing 
symmetries. This assumption is likely valid in Sr 2 Ir 04 , where most of FS is of j = 1 /2 character 
and the nearest-neighbour pseudo-spin dependent terms can be absorbed into the definition of 
the j = 1/2 states . Given the similarity in lattice structure and Mott physics between Sr 2 Ir0 4 
and La 2 Cu0 4 , it was proposed that a pseudo-spin-singlet d- wave high temperature SC phase like 
cuprates may emerge in electron-doped iridates. As discussed above, the hole-doped case may 
be significantly different due to the competition between Hund’s coupling and SOC. With finite 
Hund’s coupling one may expect an effective interaction through the exchange of ferromagnetic 
pseudo-spin fluctuations. So long as the Hund’s coupling is not so large as to drive the system into a 
ferromagnetic state, one then would expect to generate a pseudo-spin triplet pairing. Indeed, recent 
large-scale dynamical mean-field theory calculations support such a picture. In these simulations 
a rf-wave pseudo-spin-singlet SC was found on the electron-doped side, while a topological p + ip- 
wave pseudo-spin-triplet SC was found on the hole doped side . These may be relevant for hole- 
doped Sr 2 Ir0 4 , or possibly even Sr 2 Ru0 4 . A schematic phase diagram as a function of doping and 
Hund’s coupling is shown in Fig. 5b. 

Experimental studies on physics of doping holes or electrons into Sr 2 Ir0 4 have been carried 
out through a variety of means; these include electron doping through depleting oxygen 4 ”, La- 
substitution i0,51 , or surface doping 4 as well as hole-doping through substitution of Rh for Ir 50,53 . 
While SC has not been discovered, in each case modest doping suppresses the magnetic and 
strongly affects transport. Doping the bilayer Sr 3 Ir 2 07 with 5% La induces a robust metallic state, 
with the resistivity showing a rapid drop below 20 K. The magnetic order remains finite despite 
strong suppression of T c 54 . This is in contrast to the case of La-doped Sr 2 Ir0 4 , where the magnetic 
order is completely suppressed in the metallic state. In addition to the transport and magnetic 
properties, angle resolved photoemission spectroscopy (ARPES) provides interesting insights into 
physics of the doped iridates. It was reported that Rh-doped Sr 2 Ir0 4 exhibits a pseudo gap and 
Fermi arcs 4 similar to the doped cuprates. Similar Fermi arcs were also found in surface electron 
doped Sr 2 Ir0 4 along with the anomalous waterfall-like feature in the ARPES spectra . Though 
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similar Fermi arc-like features observed in La-doped Sr 3 Ir 2 07 , have been attributed to Fermi pock¬ 
ets masked by matrix element effects 6 , this does not explain the temperature dependence of the 
Fermi arcs in Ref. [52], suggesting a different origin. 


IV. SPIN LIQUIDS AND UNCONVENTIONAL MAGNETIC ORDERS IN HONEYCOMB 

IRIDATES 

Despite tremendous efforts, the search for a quantum spin liquid in a real material remains 
unresolved. A significant amount of attention in this search has been directed towards geometri¬ 
cally frustrated antiferromagnets, such as those on Kagome or triangular lattices. Many of these 
systems are described by Heisenberg-like models, possibly extended with ring exchange terms or 
with small anisotropies. While some candidates exist, the combination of theoretical and experi¬ 
mental uncertainties has made definite confirmation of the spin liquid state difficult. A promising 
approach to circumvent some of these theoretical difficulties is based on Kitaev’s exactly solvable 
honeycomb modeP . This highly anisotropic compass model is frustrated not by the geometry of 
the lattice, but by the intertwining of spatial and spin degrees of freedom. One can generalize this 
approach, defining such exactly solvable models on trivalent lattices in two or three dimensions. 
The presence of the solvable point present in the model then provides a controlled starting point 
to study the stability of the spin liquid and possible nearby unusual ordered phases. The challenge 
is then to find materials that implement these models. 

Here we discuss possible realizations of these Kitaev-type models in Mott insulators with strong 
SOC. In two dimensions, there are the well-studied honeycomb iridates Na 2 Ir0 3 and ar-Li 2 Ir0 3 , 
which host magnetically ordered phases thought to be proximate to the spin liquid phase; and the 
honeycomb ruthenium chloride ar-RuCl 3 , which has been recently proposed as a Ad analogue. Two 
three-dimensional lattices have also been synthesized in polymorphs of the honeycomb or-Li 2 Ir0 3 . 
These are the hyper-honeycomb /J-Li 2 Ir0 3 and stripy-honeycomb y-Li 2 Ir0 3 . First, we start with 
an overview of their shared basic physics: edge-shared octahedra. 


A. Spin-orbit Mott insulators with edge-shared octahedra 

Finding a system that realizes Kitaev’s honeycomb model or its analogues is difficult; no sym¬ 
metry principle prevents the introduction of other interactions, such as Heisenberg exchange ’ . 
The best one can hope for is the dominance of Kitaev exchanges among symmetry allowed in¬ 
teractions. One route to achieve this was introduced in the pioneering work of Jackeli and Khal- 
iullin based on spin-orbit entangled pseudo-spins in a 90° bonding geometry. This provides a 
natural mechanism to generate bond-dependent Ising interactions , the building block of the Ki¬ 
taev model. For concreteness, we discuss this in the context of iridium oxides, but translation to 
related compounds such as a-RuCl 3 is straightforward. 

In contrast to the perovskite iridates discussed in Sec. Ill, we consider materials where the oxy- 
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gen octahedra of neighbouring iridium atoms share an edge, rather than a corner. In addition, these 
materials are farther from the itinerant limit discussed for the perovskite case, lying firmly in the 
Mott regime. We assume that the dominant exchange pathways between iridium 5 d orbitals pro¬ 
ceed through the oxygen 2 p orbitals with a large gap A dp between the 5 d and 2 p states. Integrating 
out these oxygen states yields an effective inter-orbital hopping t = t 2 dn /A dp with t dpn being the 
Slater-Koster /r-overlap and between the 5 d and 2 p orbitals. The pair of 5 d orbitals linked depends 
on the edge shared; d yz and d xz orbitals mix on z-links, d lJZ and d xy on z/-links and d xz and d xlJ on 
x-links, as shown in Fig. 6. To understand the Mott insulating phase we will consider the effective 
model in the strong coupling limit where U, Jh, T » t, which can be expressed entirely in terms of 
the j = 1/2 pseudo-spins. The simplest approach is to first project the multi-orbital model of Eq. 
(2) into isolated j - 1/2 bands prior to taking the strong-coupling limit, as was done in Sec. Ill B. 
One finds that the inter-orbital hoppings t vanish when projected into the j - 1/2 bands, and thus 
no exchange is generated. Similarly, even in the full multi-orbital problem, no exchange arises 
when Hund’s coupling is absent. These results suggest that in this case the complete multi-orbital 
nature of the problem must be considered to account for the interactions between the j - 1/2 
spins. To include both the j - 3/2 virtual states and the effects of Hund’s coupling, one must use 
a more realistic limit where U, Jh » T or U, A » J H , obtaining an effective interaction of Kitaev 
type : 


8 J ft 
3 U 2 


J y J 7 
J 1 u 2 ’ 


( 10 ) 


where y indicates the type of edge shared by the neighbouring j — 1/2 spins Ji and J 2 . 

This Hamiltonian is only applicable when the oxygen mediated processes are dominant and A dp 
is very large. If other hopping interactions such as direct 5d-5d overlap or other super-exchange 
pathways 60,61 are included, then one expects all symmetry allowed terms to be generated. For a 
pair of ideal edge-shared octahedra, the bond symmetry allows two additional terms: a Heisenberg 
interaction J and a symmetric off-diagonal exchange T 62 . An effective exchange model for the pair 
of spins is then 

yji • j 2 + kj\ji + r [ji / 2 + j\ j "), (ii) 


where a, /? indicate the two spin directions not equal to the edge-type y. If we assume the exchange 
physics is local to the pair of octahedra, then fixing the parameters on one edge is sufficient to 
determine the others related by symmetry. Details of the dependence of J, K and Y on the atomic 
interactions for two of the limiting schemes described here can be found in Rau el al. . 


B. Distortions 

The crystal structure of the honeycomb compounds Na 2 Ir0 2 and G'-LTIrO,, and the three- 
dimensional analogues /i-LFIrCF and y-Li 2 Ir 03 deviate from this idealized picture due to the 
presence of trigonal and monoclinic distortions. Following the discussion in Section IIIC, for 
trigonal distortion, in Eq. (6) one has n = (x + y + z)/ V3 or one of its equivalents. In the absence 
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of SOC, this splits the t 2g levels into an a\ g singlet and e g doublet. When SOC is present, the 
highest-lying Kramers doublet connects to the j = 1/2 states and can be found explicitly in Ref. 
[2]. Monoclinic distortion cannot lift the degeneracy any further, but will complicate the structure 
and analysis of these ground state wave-functions. Corrections to interactions between the j = 1/2 
spins will be introduced through these changes in the atomic wave-functions, as well as through 
the kinetic parts of the exchange processes which depend on the geometry of the oxygen ions. 

Rather than attempt to connect the details of the distortions to the exchange interactions directly, 
we consider the new interactions allowed by these reductions of symmetry. In the materials under 
consideration, there are two or more symmetry inequivalent sets of nearest-neighbour bonds, one 
with higher symmetry and the rest with lower symmetry. The higher bond symmetry groups 
are 222 for /TLYIrO, and y-Li 2 Ir 03 and 2/m for o-LFIrCF and Na 2 Ir 03 . For a bond with 2/m 
symmetry, the allowed exchanges include those in Eq. (11) but with an additional symmetric 
off-diagonal exchange 

r [j/jI + f\j\ + j\r 2 + /j7?). ( 1 2) 

As an example, trigonal distortion of the oxygen octahedra can introduce such a term . For a 
bond with 222 symmetry, such a term is forbidden, but due to the lack of inversion about the bond 
center an analogous DM interaction is allowed 

d (jij\ + yf j\ -J\J 2 - J]J 2 ). (13) 

The lower symmetry bonds of or-Li 2 Ir0 3 , /3-Li 2 Ir0 3 , and Na 2 Ir0 3 have only inversion about the 
bond center while that of y-Li 2 Ir0 3 has no symmetry at all. In the compounds with an inversion, 
this implies the absence of DM interaction on these bonds. 

Further neighbour exchange interactions can be analyzed in a similar fashion, though the num¬ 
ber of allowed interactions grows quickly. These are thought to be important “ in Na 2 Ir 03 where 
the Mott gap is not too large . Some discussion of the symmetry allowed interactions in Na 2 Ir 03 
for second and third nearest neighbours can be found in Yamaji et al. and Sizyuk et al. . 
While there has been some work on further neighbour interactions in the three-dimensional 
compounds 69,70 , their effects and importance remain largely unaddressed. 


C. Kitaev physics in two dimensions 

With this framework in hand, we first turn to possible realizations of Kitaev’s two-dimensional 
honeycomb model. Promising candidates thought to be proximate to this physics are the honey¬ 
comb iridates Na 2 Ir 03 and a-LFIrC^, and more recently, o'-RuC^. These materials do not have a 
spin liquid ground state, but instead order magnetically as temperature is lowered. To understand if 
the proximity to Kitaev physics is governing their behaviour we first study these magnetic phases. 
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(a) Honeycomb plane 


(b) Unit cell of Na 2 lr 03 


FIG. 6: Crystal structure of the layered Na 2 Ir0 3 . (a) The unit cell of Na 2 Ir 03 showing the 
quasi-two dimensional layered structure, (b) The honeycomb plane of Na 2 Ir 03 showing the 
labeling of the nearest neighbour bonds, with a y = x, y, z type bond labeled as a/3(y). 


1. Experimental review 

The crystal structures of these compounds (see Fig. 6) takes the form of layers of edge-shared 
Ir06 octahedra arranged in a honeycomb lattice 1, . These octahedra are compressed in the direc¬ 
tion perpendicular to the honeycomb planes, with further monoclinic distortions along a preferred 
axis ’ . A 340 meV optical gap and estimates of 3T/2 ~ 0.5 eV for SOC identify Na 2 Ir0 3 
as a candidate for a j = 1 /2 spin-orbit Mott insulator. Fits of the magnetic susceptibility confirm 
the effective spin-1/2 picture, giving a magnetic moment ~ 1.82 p B with large antiferromagnetic 
Curie-Weiss temperature 0 C w-116 K . The low temperature antiferromagnetic ordering tran¬ 

sition seen near T N ~ 15 K suggests substantial frustration. Neutron scattering and resonant X-ray 
scattering studies ■ ■ have unambiguously identified this ordering as having a zigzag structure. 
This zigzag order consists of spins aligned ferromagnetically along one direction, with the chains 
alternating antiferromagnetically as illustrated in Fig. 7. Dynamical probes, such as inelastic 
neutron scattering (INS) and resonant inelastic X-ray scattering (RIXS), have provided additional 
clues via access to their excitations. RIXS studies have shown that there is a branch of magnetic 
excitations at high-energy near ~ 30 meV. Though some dispersion can be identified, any infor¬ 
mation about the excitations below ~ 10 meV is lost due to limitations in energy resolution. INS 
studies offer a complementary picture, but due to the large neutron absorption cross-section of Ir, 
they have been limited to powder samples 7 providing a higher energy resolution but only a few 
details of the distribution in wave-vector. Nonetheless, below the Neel temperature, two key fea¬ 
tures can be resolved: scattering near the magnitude of the zigzag ordering wave-vector is present 
down to at least 2 meV and there is an absence of scattering at small wave-vectors and energy. Re- 
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cent diffuse magnetic X-ray scattering in the paramagnetic phase has provided an experimental 
confirmation of dominant Kitaev interactions, validating the theoretical arguments of Sec. IV A. 
These measurements show a clear locking of the spatial and spin orientations characteristic of the 
bond-dependent Kitaev exchange. For ar-Fi 2 Ir0 3 the situation is less clear. Experimentally, an an¬ 
tiferromagnetic ordering transition is seen at ~ 15 K, as in Na 2 Ir0 3 , but a- Fi 2 Ir0 3 has a smaller 

Curie-Weiss temperature of 0 C w-33 K. While there have been reports of an incommensurate 

ordering wave vector lying in the first Brillouin zone 79 , the details of the magnetic ordering pattern 
and the structure of the low energy excitations remain largely unresolved. 

A number of studies have tried to elucidate the properties of these materials more indirectly 
through elemental substitution. One promising approach is to dope isoelectronically from Na 2 Ir0 3 
to Q'-Li 2 Ir0 3 as (Na 1 . x Li x ) 2 Ir0 3 . For x < 0.25, uniform solid solutions can be obtained, with both 
T N and |(9 CW | suppressed with increased doping ■ . While a quantum critical point near x c ~ 0.75 
was found in Ref. [80], indications of phase separation reported for the range 0.25 < jc < 0.6 
complicate this identification. Further study is needed for x > 0.6 near the ar-Fi 2 Ir0 3 end of this 
range to clarify the issue. Another approach is dilution of Ir 4+ with non-magnetic Ti 4+ , forming 
Na 2 (Iri_ x Ti x )0 3 or Fi 2 (Iri_ x Ti x )0 3 . As the magnetic lattice of Ir is depleted, both systems enter a 
spin-glass phase, consistent with the high frustration indicated by T n /\6qw\- In both compounds 
the spin glass ordering temperature T g decreases roughly linearly as a function of x until to the 
site-percolation threshold at x p ~ 0.3. These two cases are distinguished by the dependence of 
the Curie-Weiss temperature on dilution: the Na 2 (Iri_ x Ti x )0 3 shows a marked decrease as x is 
increased while for the Fi 2 (Iri_ x Ti x )0 3 it is essentially constant. This has been interpreted as 
evidence of more significant long range interactions in a-Fi 2 Ir0 3 compared to Na 2 Ir0 3 . 


2. Theory 

The balance of experimental evidence has suggested that large Kitaev interactions are necessary 
to understand the physics of Na 2 Ir0 3 and a-Fi 2 Ir0 3 . However, the low-energy physics and the 
ground state selection are dependent on the details of the perturbations that take us away from 
the Kitaev limit. These questions have led to a number of theoretical proposals that have been 
put forth to explain the appearance of zigzag ordering in Na 2 Ir0 3 and the nature of the magnetic 
ordering in ar-Fi 2 Ir0 3 . 

Early work on Na 2 Ir0 3 primarily focused on the Heisenberg-Kitaev (HK) model and its ex¬ 
tensions. In the HK model, the Kitaev exchange is supplemented by a conventional Heisenberg 
interaction 

<i4) 

This model has since attracted much theoretical attention, not only due to its relation to the hon¬ 
eycomb iridates, but also as a model system to study the stability of the Kitaev spin liquid and its 
neighboring phases 8 ' . The remainder of the phase diagram of this model can be understood al¬ 

most completely with the help of a four-sublattice spin rotation—the so-called Klein duality. This 
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FIG. 7: Phase diagram of the Heisenberg-Kitaev model. The model is parametrized as J = cos (p, 
K = sin (p , as discussed in the main text. Squares and diamonds show points related by the Klein 
duality transformation. Phases of the Heisenberg-Kitaev model are shown inside, these are the 
ferromagnet (FM), antiferromagnet (AFM), stripy (ST) and zigzag (ZZ). 


duality maps the HK model to itself, but with the modified parameters J' = - J and K' - K + 27. 
Applying this duality to the ferromagnetic and antiferromagnetic Heisenberg limits yields two new 
well-understood limits. These are the stripy (ST) phase at K = -2 J < 0 dual to the ferromagnet 
(FM) and the zigzag (ZZ) phase at K = 2J > 0 dual to the antiferromagnet (AFM). These ordered 
phases are illustrated in Fig 7. The Kitaev spin liquid is stable to finite J and maps to itself under 
the Klein duality. These five phases encompass all phases present in the model. This has been 
borne out using a number of theoretical methods 60,8 ,M ; the full phase diagram parametrized as 
J = cos (p, K = sin0 is presented in Fig. 7. 

For the presumed ferromagnetic Kitaev interaction (K < 0), the neighboring magnetic phases 
are a ferromagnet and a stripy phase 6 . It is possible to obtain a zigzag ground state with this 
model, but only when the Kitaev interaction is antiferromagnetic (K > 0) and the Heisenberg 
interaction is ferromagnetic (J < 0). The need for a large antiferromagnetic Kitaev interaction 
is at odds with expectations from a picture of dominant oxygen-mediated super-exchange. Some 
efforts have been made to justify this parameter regime, while others 6 ’ ■ have added addi¬ 
tional couplings to the ferromagnetic Kitaev limit in an attempt to resolve the discrepancy. Recent 
ab-initio calculations 65, have supported the view that ferromagnetic Kitaev interactions are the 
proper starting point. Starting from this Kitaev limit, theoretical studies have considered the ef- 
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fects of the symmetry allowed nearest-neighbour interactions - such as the T from Eq. (11) 
and F from Eq. (12), as well as further neighbour interactions ' . While the details differ, in all 
of the above works a zigzag phase can be stabilized near the ferromagnetic Kitaev limit by some 
combination of further neighbour and anisotropic interactions. 

Even with the scarcity of experimental information, several theories have been put forth to 
explain the possibility of incommensurate ordering in a'-LEIrCE. A common approach is to move 
away from some of the limiting regimes studied for Na 2 Ir 03 . For example, in Ref. [84], the 
HK model was extended to include second nearest-neighbour Heisenberg and Kitaev exchange 
interactions. In this regime one can find an incommensurate spiral phase with wave-vector lying 
in the first Brillouin zone, as has been reported in a- LkMk. In Refs. [63] and [88], one perturbs 
away from the zigzag phase stabilized by anisotropic exchanges Y and P, which is relevant for 
Na 2 Ir0 3 , to find similar incommensurate phases. Other approaches begin from a more monoclinic 
limit, considering quasi-one-dimensional chains or strengthening Kitaev interactions on one set 
of bonds 6 . The lack of details on the reported incommensurate spiral phase coupled with the same 
large parameter space that plagues Na 2 Ir0 3 leaves many questions open for o'-Li 2 Ir0 3 . 


3. a-RuCl3 

The 4 d compound ar-RuCk has recently attracted attention as another possible system to ex¬ 
plore Kitaev physics. As in honeycomb iridates, the Ru 3+ ion has a d 5 configuration but with an 
octahedral cage of Cl“ rather than O 2 . These octahedra are then arranged in a layered, edge- 
shared honeycomb network with space group P3 1 12 . While first identified as a band insulator 91 , 
spectroscopic studies 9 ■ have since converged to a picture of ar-RuCk as a Mott insulator. In¬ 
deed, this is supported by measurements of the magnetic susceptibility, which takes a Curie-Weiss 
form with an average moment ~ 2.0 - 2.3 not too far from what would be expected for a 
spin-1/2 ■ . Going to lower temperatures, antiferromagnetic order appears near ~ 7 K, as seen in 
specific heat and magnetic susceptibility measurements 1-1 . Neutron scattering experiments ■ 
have identified this order as being zigzag 7,99 as seen in Na 2 Ir0 3 . Further, a broad specific heat 
feature seen near ~ 15 K • may be an indication of a multi-step transition process. 

While strongly suggestive of the same physics seen in the honeycomb iridates, one may be 
concerned with the applicability of the j = 112 picture itself. Given the considerably smaller 
SOC in the 4 d orbitals of Ru 3+10 °, strong overlap of the j = 1/2 and j = 3/2 bands is possible 
and would require a more complicated description of the Mott insulating phase. One proposed 
resolution , motivated by ab-initio calculations, invokes the enhancement of SOC by electronic 
correlation effects. This SOC enhancement cleanly separates the j = 1/2 states from the j = 3/2 
and increases the tendency towards a Mott phase. Several theoretical descriptions ' 6,99,11 have been 
put forth to explain the appropriate low-energy physics, as well as the nature of the ordered phase 
that appears at low temperature. 

Some of the key differences between a-RuCk and the honeycomb iridates may render it a 
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(a) /?-Li 2 Ir0 3 


(b) y-Li 2 Ir0 3 


FIG. 8: Crystal structures of the (a) hyper-honeycomb /J-Li 2 Ir0 3 and the (b) stripy-honeycomb 
y-Li 2 Ir0 3 . Each structure is formed from a network of approximately edge-shared Ir0 6 
octahedra, with bonds of type y - x, y, z labeled as nfi(y). 


much better system to study Kitaev magnetism. Indeed, geometrically the local environment of the 
Ru ion is nearly distortion-free 93, , lacking the significant trigonal or monoclinic distortions that 
have been reported ■ 2 in Na 2 Ir0 3 and a- Li 2 Ir0 3 . Access to the details of the magnetic excitation 
spectrum, due to the absence of neutron absorbing Ir atoms 9 , should further illuminate the physics 
of this material. Aside from pinning down the theoretical description, with more information and 
fewer parameters it may be easier to perturb the system toward the sought-after Kitaev spin liquid 
phase . 


D. Kitaev physics in three dimensions 

The search for spin liquids in three dimensional materials has focused on a handful of frus¬ 
trated lattices, for example on the pyrochlore ! ■ or the hyper-kagome, realized in the iridate 
Na 4 Ir 3 0 8 ■ . As in the two-dimensional cases discussed in Section IV have been put forth, 

the definite confirmation of the spin liquid state remains very challenging both theoretically and 
experimentally. While further work is warranted, following our discussion in two-dimensions, an 
alternative route to a three-dimensional spin liquid might lie through an extension of the Kitaev’s 
honeycomb model. 

The exact solvability of Kitaev’s 2D honeycomb model depends on the three-fold coordination 
of each site and the appropriate assignment of Ising components among the bonds, i.e. the Ising 
interactions are along the x, y, or z directions while bonds with the same Ising components do not 
share the same site. These properties can be realized on many lattices, hence Kitaev’s 2D honey¬ 
comb model can be generalized to three-dimensions. One such example is a spin-1/2 model de¬ 
fined on a 3D deleted-cubic lattice, whose exact spin liquid solution and excitations were studied in 
Ref. [106]. Fortuitously, an iridate material with the Ir ions residing on a topologically-equivalent 
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lattice was synthesized several years later: the hyper-honeycomb /TLMrCV. Almost simultane¬ 

ously, another polymorph—the stripy-honeycomb y-LblrO^—was discovered independently. 
Both polymorphs contain the correct trivalent lattice structure, bond geometry, and large SOC 
to possibly realize a Kitaev model, thereby providing a new avenue toward the discovery of a 
3D spin liquid. However, magnetic orders were experimentally observed in both /i-LFIrCF and 
y-Li 2 Ir 03. 10 ’’ Similar to the two-dimensional case, both theoretical and experimental work has 
focused on the magnetic phases in order to understand the role played by Kitaev physics. Addition¬ 
ally, theoretical investigations have been undertaken to explore this new 3D spin liquid appearing 
in the Kitaev limit. 


1. Experimental review 

The hyper- and stripy-honeycomb iridates are built from networks of edge-shared octahedra 
much like the 2D honeycomb iridates, but these networks are intrinsically three-dimensional. Both 
of these materials are insulating, with the magnetic susceptibility indicating a moment size close 
to that of a spin-1/2 degree of freedom. As in the case of Na 2 Ir 03 and o'-LMrCF, this is con¬ 
sistent with a picture of localized j = 112 moments. An antiferromagnetic ordering transition 
is seen at ~ 38 K in both materials through a kink in the magnetic susceptibility and a peak in 
the specific heat, though Curie-Weiss fits demonstrate predominant ferromagnetic exchanges in 
both compounds. In the case of y-Li 2 Ir 03 , these exchanges are believed to be anisotropic 
as torque magnetometry measurements have revealed a temperature dependence in the anisotropy 
of the magnetic susceptibility. Further studies of these orderings using magnetic resonant x-ray 
scattering have identified the magnetic states in both the /TLi 2 Ir 03 and y-Li 2 Ir 03 as incommensu¬ 
rate, non-coplanar, counter-rotating spirals, as illustrated in Fig. 9. Particularly striking is the close 
agreement of the incommensurate propagation vectors Q ~ [0.57,0,0] in both compounds. 109,110 
While closely related, these two ordering patterns have subtle differences: In /TLFIrO,, the mo¬ 
ments projected along each of the three orthorhombic directions, a, b, and c transform under the 
same irreducible representation for a magnetic structure with wave-vector [ q , 0,0]. In contrast, the 
a and c components of the spiral phase of y-Li 2 Ir0 3 transform under a different irreducible repre¬ 
sentation from that of the b component. The similarity of the orderings, ordering temperature, as 
well as local geometry of these crystals suggest that a common understanding might be possible. 


2. Theory 

These experimental discoveries have motivated a number of theoretical proposals to explain 
and explore these 3D honeycomb materials. ’ ~ Several works attack /TLFIrCF and y-Li 2 Ir 03 
directly, starting from a ferromagnetic Kitaev limit and considering the effects of further perturba¬ 
tions. These include the allowed nearest-neighbour exchanges discussed in Section IV A, as well 
as those between further neighbours and those induced by the monoclinic nature of the lattice. 
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(b) Ordering in y-LFIrOs 



FIG. 9: Illustration of the incommensurate, non-coplanar, counter-rotating spirals in (a) 
jS-LFIrCV' and (b) y-LFIrOi 11 . For both these materials, the ordering wave-vector is aligned to 
the orthorhombic a direction and three unit cells along that direction are depicted. The 
non-coplanar nature is readily seen in the perspective illustrations on the left while the 
counter-rotating nature of these spiral orders are evident in the projected illustration on the right. 
Color indicates whether the b component of the spin is positive (red), negative (blue) or zero 
(black). Due to the difference in crystal symmetry of the two lattices, the b component of the 
moments in y-Li 2 Ir 03 transforms differently from the a and c components. 


Using these interactions, several routes have been identified to stabilize an incommensurate spiral 
ordering near the ferromagnetic Kitaev limit ■ ■ . Though capable of stabilizing the ground 

state, these theories have not been tested against current experimental findings such as magneti¬ 
zation, torque magnetometry, and thermodynamic measurements. Moreover, experiments such as 
inelastic scattering that probe low-energy dynamics of these systems have yet to be conducted. 
These comparisons and results should provide a more comprehensive understanding of the role of 
Kitaev physics in the unconventional magnetism encountered in these 3D honeycomb iridates. 
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In addition to work related to the experimental magnetic phase, theoretical work has also 
focused on the spin liquids that would be relevant if the Kitaev limit is stabilized on these 
lattices. 6,111,112,11 - Identical to the 2D honeycomb case, the exact solution takes the form 
of free Majorana fermions in the presence of a static background Z 2 gauge field, and both gapped 
and gapless phases can be accessed by tuning the bond anisotropy away from the isotropic limit. 
Unlike the 2D case, the lack of mirror symmetry implies that the background flux generated by 
the gauge field in the ground state may not be uniform. Further, due to the additional constraints 
of the three-dimensional geometry, flux excitations of the Z 2 gauge field are restricted to take 
the form of loops. These loops play a key role in the finite temperature properties of the model, 
driving a phase transition that separates the low temperature spin liquid and high temperature 
paramagnet . Near the isotropic limit where the phase is gapless, the spin liquid possesses un¬ 
usual nodal lines (co-dimension 2) of gapless Majorana modes. These nodal lines are topologically 
stable and, due to the bulk-boundary correspondence, induce gapless modes on the surface. 1 In 
the presence of time-reversal symmetry breaking, these line nodes become topologically protected 
Weyl points with associated surface Fermi arcs. 1 

Aside from the hyper- and stripy-honeycomb lattices, a number of other hypothetical triva- 
lent lattices in three-dimensions have also been considered. The family of lattices termed the 
harmonic-honeycombs is a generalization of the hyper- and stripy-honeycomb lattices. These 
lattices are constructed by exploiting the fact that the hyper- and stripy-honeycombs can be con¬ 
sidered as rows of honeycomb lattices stacked in an alternating fashion. By varying the number 
of complete honeycomb rows in such a construction, the harmonic-honeycomb series of lattices 
is generated. Each lattice of the series is denoted as H (N), where N is an integer that refers 
to the number of complete honeycomb rows used in the construction. The hyper- and stripy- 
honeycomb lattices are ^(0) and H ( 1) respectively in this notation. Like the hyper- and stripy- 
honeycomb lattices, the Kitaev spin liquids in all finite-A harmonic-honeycomb lattices possess 
topologically protected line nodes of gapless Majorana fermions and surface gapless modes. 1 In 
contrast, another 3D trivalent lattice—the hyper-octagon lattice—possesses a gapless phase with a 
2D Majorana Fermi surface. All of these theoretical lattices can be embedded in structures con¬ 
taining edge-shared Ir0 6 octahedra, hinting at the possibility of realizing a 3D Kitaev spin liquid 
in a yet-to-be-discovered material. 


V. OUTLOOK 

The wide variety of phenomena observed and predicted in iridates and related materials has 
generated immense activity in the field. Our objective in this review was not to be exhaustive but 
rather to instill a sense of breadth by highlighting a few specific directions recent studies have 
pursued. In particular, we examined the Ruddlesden-Popper series of iridates and their related 
heterostructures where unconventional magnetism, topological band structures, and superconduc¬ 
tivity have been seen or predicted. We also examined honeycomb materials, both in two and three 
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dimensions, where the joined effects of SOC, electron correlation, and lattice geometry give rise 
to Kitaev physics and its associated magnetism. These topics merely scratch the surface of possi¬ 
bilities in the study of the interplay of SOC and electronic correlations; related research directions, 
both theoretical and experimental, further illustrate the richness of this field. With this rapid de¬ 
velopment, however, there remain several obstacles that if resolved, can substantially enhance our 
understanding of SOC effects and propel the community forward. 

One issue that has hindered the study of iridate materials is the large neutron absorption cross- 
section of Ir. This property has restricted neutron scattering experiments to powder samples or 
large single crystals, where the former studies lack wave-vector details while the latter studies 
may require years of growth technique refinement. Although great headway has been achieved 
using resonant X-ray scattering as an alternative probe, higher energy resolution and the ability to 
perturb the measured system in situ would enable these techniques to be a true alternative to neu¬ 
tron techniques. Sample quality is also in need of continual refinement in many of the considered 
compounds, especially the honeycomb iridates where Na and Li readily oxidizes, as discussed in 
Sec. IV C 1. ffomogeneous solutions of single crystal of (Nai. x Li x ) 2 Ir0 3 over a wide range of x 
(and especially x = 1) would greatly enhance our understanding of Kitaev physics and its con¬ 
nection with microscopic details of materials. In addition to elemental substitution, other methods 
to perturb these systems toward the Kitaev regime or MIT, for example with applied pressure, 
could also provide ample insight. Even with the well-studied perovskites iridates, there are still 
opportunities for major breakthroughs. Direct evidence of superconductivity in Sr 2 Ir 04 or other 
iridium oxides is still intensely pursued, as is evidence and examples of topological phases. Aside 
from strontium iridates, other avenues to realize topological phases have also been sought. As 
mentioned in Sec. HID, engineered phases like superlattices or thin films of 4 d or 5 d TMOs 
play a special role in this area of exploration due to their tunability of SOC, correlation, and 
lattice distortions. In this regard, materials like thin film pyrochlore iridates 1 ~ and epitaxial 
perovskites ' ’ _ have received considerable attention. On the theoretical front, extending our 
understanding of topological phases of Sec. Ill D beyond the non-interacting limit is an important 
goal. A large effort has focused on the consequences of interactions and disorder in Weyl semi¬ 
metals, especially in relation to its characterization , transport properties 128,129 , and instabilities 
towards exotic phases . Related is the study of Coulomb interaction on the quadratic band 
touching in pyrochlore iridates, which can induce a quantum critical non-Fermi liquid, a Weyl 
semi-metal, or a topological insulator depending on the preserved symmetries. These are but 
some examples where theoretical predictions of interacting topological phases may find realiza¬ 
tions in iridate materials, and continual search for these systems will be a central focus in the years 
coming. 

In addition to addressing outstanding issues, the discovery of new compounds has also driven 
promising lines of investigation. These new materials provide a wonderful arena both to test theo¬ 
retical predictions and to discover new phenomena related to strong SOC. A recent example is the 
Mott-insulating Ba 3 IrTi 2 0 9 that does not order magnetically. It has been argued to realize the 
ffK model on a triangular lattice, leading to frustration in both geometry and exchanges, which 
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can lead to a variety of exotic phases. ’ Another recently synthesized compound thought to host 
the HK model is the honeycomb LBRI 1 O 3 , which is closely related to the 2D honeycomb materials 
discussed in Sec. IV C. This 4 d 5 material does not magnetically order but shows spin-glass be¬ 
havior that may be related to disorder or stacking faults. With this reported rhodate, alluring ideas 
like the synthesis of 3D honeycomb rhodates or isoelectronic substitution of honeycomb materials 
will allow the tuning between 4 d and 5 d materials to tilt the balance between SOC, crystal field, 
and correlation effects. 

Intimately related to lattice geometry is how the Ir06 octahedra are inter-connected; the ma¬ 
nipulation of this connectivity has also been an interesting approach to generate new compounds. 
Our review has focused on comer- and edge-shared octahedra, but there are a variety of com¬ 
pounds with face-sharing octahedra where SOC plays an important role in the behavior of the 
material. The Mott insulating and magnetically ordered double perovskites La 2 ZnIr06 and 
La 2 MgIrO A serve as yet another example 1 ' : the Ir0 6 octahedra in these double perovskites are 
spatially separated because they do not share corners, edges, nor faces. The resulting reduction 
in orbital overlaps between nearest neighbor 5 d Ir orbitals leads to a decrease in bandwidth or, 
equivalently, a relative enhancement of SOC and correlation effects. By further spatially isolating 
the Ir0 6 octahedra, the validity of the localized picture of j = 1 /2 orbitals can be probed. The 
magnetically ordered compound Ca 4 Ir0 6 is one such example . 1 ' The large separation between 
octahedra in this compound and near-ideal Ir 06 geometry has led to one of the first detections of 
distortion-free j = 1/2 orbitals in this compound. The compound Sr^Lan_ v Ir 4 0 2 4 has also been 
used to test the j = 1/2 states but offers an additional twist: with the tuning of 1 < x < 5, the 
valence of the iridium ions can be changed between Ir 4+ ( d 5 ) and Ir 5+ (J 4 ). In the strong SOC 
limit, the d 4 filling yields a non-magnetic state since the atomic j - 3/2 states are completely oc¬ 
cupied and the j - 1/2 orbitals remain unfilled. Measurements of the magnetic properties of this 
series of materials indeed show a transition between the j = 1/2 and non-magnetic state, indicating 
that strong SOC is essential in the description of spatially isolated Ir0 6 octahedra. However, this 
strong SOC limit may not apply to other compounds with the d 4 configuration: examples include 
the layered perovskite Ca 2 Ru0 4 , the honeycomb materials /OR 11 O 3 (A=Li, Na) 143 , and yet an¬ 
other variety of double perovskites, Sr 2 YIr0 6 and La 2 MRu0 6 (M=Mg and Zn ) 14 . Magnetic 
transitions have been observed in these compounds and the origin of such magnetism is under ac¬ 
tive investigation. Proposals such as Van Vleck-type Mott insulators 146 , strong non-cubic crystal 
field effects 14 , and competition between super-exchange and SOC have been investigated in the 
context of these magnetic d 4 systems. 

Beyond integer fillings, systems with partial filling have also received attention and the multi¬ 
orbital nature of these transition metal compounds may play an elevated role. In the recently 
synthesized Ba 5 AlIr 2 0n, dimers of face-shared Ir06 octahedra have an average valence of Ir 4 5+ , 
yet the material is a Mott insulator . 1 In the thiospinel CuIr 2 S 4 , where sulphur plays the role of 
oxygen in the octahedral cages, the average valence is Ir 3 - 5+ .' Several recent studies have shown 
that SOC, which have been neglected in earlier works, may be important in the newly observed 
low temperature paramagnetic state. Its cousin, CuIr 2 Se 4 , is a metal, but when doped with Pt, 
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shows evidence of superconductivity: an exciting opportunity considering the local structural 
similarities between IrSe 6 and the Ir06 octahedra discussed throughout this article. Another ex¬ 
ample is the hyper-kagome, Na 3 Ir 3 0 8 , which has an average valence of Ir 433+ and thus can be 
considered as a 1/3-doped version of the spin-liquid candidate Na 4 Ir 3 0 8 mentioned briefly in Sec. 
IVD. 1 Unlike the Mott-insulating Na 4 Ir 3 0 8 , this new material is semi-metallic, which was ar¬ 
gued to stem from the presence of strong SOC and distortion-induced molecular orbitals. 1 With 
a combination of its chiral crystal structure and frustrated lattice, it was suggested that this semi- 
metallic material may harbour a unique platform to study non-trivial electronic transport due to 
topological effects arising from strong SOC. 

As we have seen in these brief examples, the study of the combined effects of strong SOC and 
correlations has grown to encompass a wide range of materials and phenomena. This breadth is an 
indication of the rapid development of this field, which boosts the prospects of discovering novel 
and exciting physics. With the eventual advancement in experimental, numerical, and theoretical 
techniques, many obstacles currently faced can be tackled, and perhaps even resolved, thus paving 
way for new and exciting research directions in the near future. 
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